Josephson Effect due to Odd-frequency Pairs in Diffusive Half Metals 
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The Josephson effect in superconductor / diffusive ferromagnet / superconductor (SFS) junctions is studied 
using the recursive Green function method in the regime of large exchange energy in a ferromagnet. Motivated 
by recent experiment [R. S. Keizer, et. al. Nature 439, 825 (2006)] we also address the case of superconductor 
/ diffusive half metal / superconductor junctions. The pairing function in spin-singlet and triplet channels, 
the Josephson current and their mesoscopic fluctuations are calculated. We show that the spin-flip scattering 
at the junction interfaces opens the Josephson channel of the odd-frequency spin-triplet Cooper pairs. As a 
consequence, the local density of states in half metals has a large peak at the Fermi energy. Therefore odd- 
frequency pairs can be detected experimentally by using the scanning tunneling spectroscopy. 

PACS numbers: 74.50.+r, 74.25.Fy,74.70.Tx 



Ferromagnetism and spin-singlet superconductivity are 
competing orders against each other because the exchange 
field breaks down the spin-singlet pairs. The Cooper pairs, 
however, do not always disappear under exchange fields. 
The Fulde-Ferrell-Larkin-Ovchinniko\ii2, (FFLO) state and 
the proximity effect in ferromagnets^ii^f^iiSiSiiS are typical ex- 
amples. Under exchange fields, the pairing function oscillates 
and changes its sign in the real space. As a consequence, 
superconductor / ferromagnet / superconductor (SFS) junc- 
tions undergo the O-tt transition with varying length of a fer- 
romagnet or temperature. It was also predicted that the odd- 
frequency spin-triplet pairs appear in weakly polarized ferro- 
magnets with rotating magnetization direction near the junc- 
tion interface^. These pairs have long range penetrarion into 
ferromagnets^'^. Thus Cooper pairs change their original face 
to survive under exchange fields. 

Half metal is an extreme case of completely spin polarized 
material because its electronic structure is insulating for one 
spin direction and metallic for the other At a simple thought, 
the spin-singlet Cooper pairs would not be able to penetrate 
into half metals. However recent experiment" showed the 
existence of the Josephson coupling in superconductor / half 
metal / superconductor (S/HM/S) junctions, where NbTiN 
was used as a s-wave superconducting electrode and Cr02 
as a half metal. Thus one has to seek a new state of Cooper 
pairs in half metals attached to spin-singlet superconductors. 
Prior to the experiment", Eschrig et. al}^ have addressed 
this challenging issue. In the clean limit, they showed that 
the p-wave spin-triplet pairs induced by the spin-flip scatter- 
ing at the interface can carry the Josephson current. In real 
S/HM/S junctions, however, half metals are in the diffusive 
transport regime; the elastic mean free path of an electron is 
much smaller than the size of a half metal. In addition, real 
S/HM/S junctions are close to the dirty limit because the co- 
herence length in a ferromagnet may become comparable to 
the mean free path. In such diffusive half metals, the p-wave 
symmetry of Cooper pairs is not possible because the pair 
wave function is isotropic in the momentum space due to im- 
purity scattering'^. In this paper, we provide a general theory 



for the Josephson effect in diffusive SFS junctions with arbi- 
trary magnitude of the exchange field Vex- When Vex is much 
larger than the superconducting pair potential at zero tempera- 
ture Aq, mesoscopic fluctuations of the Josephson current are 
much larger than the ensemble averaged value. In addition, 
we focus on interesting case of diffusive S/HM/S junctions. 
We show that the odd-frequency spin-triplet s-wave pairing 
state is realized in half metals and propose an experimental 
method to detect this property. 

Let us consider the two-dimensional tight-binding model 
for a SFS junction as shown in Fig. Q^a). The vector r = 
jx + my points to a lattice site, where x and y are unit vec- 
tors in the x and y directions, respectively. In the y direction, 
we apply the periodic boundary condition for the number of 
lattice sites being W . Electronic states in superconducting 
junctions are described by the mean-field Hamiltonian 
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hj. j.1 = \_^t6^r-r'\,l + i^r ^ M + 4i)(5.r,r'] O'O 

-Vir)-&Sr,r', (2) 

with Cr — (cr,t, Cr,i), whcrc cl ^ (c^ ^) is the creation (an- 
nihilation) operator of an electron at r with spin cr =( t or 
i ), c means the transpose of c, ai for ^ = 1 — 3 are the 
Pauli's matrices, and cto is 2 x 2 unit matrix. The hopping 
integral t is considered among nearest neighbor sites in both 
superconductors and ferromagnets. In a ferromagnet, the on- 
site scattering potentials are given randomly in the range of 
— V//2 < f r < Vj/2 and the uniform exchange potential is 
given hy V{r) — Vex^s, where e/ for Z = 1 — 3 is unit vector 
in a spin space. The Fermi energy is set to be 2t in a normal 
metal with Vex = 0, while a ferromagnet and a half metal are 
respectively described by Vex/t = 1 and 2.5 in Fig.Qb). The 
spin-flip scatterings are introduced at j = 1,2, — 1, and 
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FIG. 1: (Color online) (a): A schematic figure of a SFS junction 
on the tight-binding lattice, (b): The density of states for each spin 
direction. The Josephson junction is of the SNS, SFS, and S/HM/S 
type for Vex/t = 0, 1 and 2.5, respectively. 



Ltv, where we choose V{r) = Vse2- In superconductors we 
take Cr — Q and choose A = iA(T2, where A is the amplitude 
of the pair potential in the s-wave symmetry channel. 

The Hamiltonian is diagonalized by the Bogoliubov trans- 
formation and the Bogoliubov-de Gennes (BdG) equa- 
tion is numerically solved by the recursive Green function 
method''* '^. We calculate the Matsubara Green function, 

r (rr']-f /--.(^'^') ^ 

where ujn ~ (2rt + I jirT is the Matsubara frequency, n is 
an integer number, and T is a temperature. The Josephson 
current is given by 

w 

J = -ietT ^ Tr [G^„ {r', r) - (r , r')] (4) 

uJn m— 1 

with r' = r + X. In this paper, 2x2 and 4x4 matrices 
are indicated by • ~ • and • " •, respectively. The quasiclassical 
Green function method is a powerful tool to study the proxim- 
ity effect. However the quasiclassical Green function cannot 
be constructed in a half metal because the Fermi energy is 
no longer much larger than the pair potential for one spin di- 
rection. On the other hand, there is no such difficulty in our 
method. In addition, it is possible to obtain the ensemble av- 
erage of the Josephson current (J) = (l/Ns) X^tfi -^i '^'^'^ 
fluctuations 5J = \J ( J^) — ( J)^ after calculating Josephson 
cuiTent in a large number of samples Ng with different impu- 
rity configurations. These are the advantages of the recursive 
Green function method. Throughout this paper we fix the fol- 
lowing parameters: Ln — 74, W — 25, ^ = 2t, Vj = 2t 



and Ao = 0.005i (W]. This parameter choice corresponds to 
the diffusive transport regime in the N, F and HM layers. The 
results presented below are not sensitive to variations of these 
parameters. 

We first discuss the Josephson cuiTent in SFS junctions as 
shown in Fig. |3a) for T ~ O.lTc where Tc is the transition 
temperature. We assume that the spin-flip scattering at the in- 
terfaces is absent (i.e., Vs — 0) and fix the phase difference 
across the junctions ip equal to 7r/2. The results are normal- 
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FIG. 2: (Color online) (a) The Josephson current versus the ex- 
change potential Vex- At Vex ~ 2t ferromagnets become half- 
metallic as indicated by an arrow, (b) The pairing function of spin- 
singlet pairs and its fluctuations versus position j for Vex/t ~ 
(SNS) and 2.5 (S/HM/S); (c) for Vex/t = 1 (SFS). (d) The pairing 
function in three different samples for Vex/t — 1, the vertical axis 
is offset by 0.5 as indicated by the horizontal lines. The spin-flip 
scattering is absent in all panels, Vs = 0. 

ized by ( Jq) which is the Josephson current in superconductor 
/ normal metal / superconductor (SNS) junctions (i.e.. Vex = 
0). We define the coherence length = ^ D/2Vex/ao mea- 
sured in units of lattice constant oq with D being the diffu- 
sion coefficient. The Josephson cuiTent oscillates as a func- 
tion of Vex and changes its sign almost periodically. The sign 
changes of (J) correspond to the O-tt transitions in an SFS 
junction. At the same time, the amplitude of (J) decreases 
rapidly with increasing Vex- We should pay attention to the 
relation (J) <C SJ which means that the Josephson current is 
not the self-averaging quantity. It is impossible to predict the 
Josephson cuiTent in a single sample J,; from the ensemble 
average (J) because Ji strongly depends on a microscopic 
impurity configuration. In fact, the Josephson current flows 
in a single sample even if (J) = at the transition points. 
Roughly speaking, (J) vanishes because half of samples are 
0-junctions and the rest are Tr-junctions'^. Since (J) = 0, 
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6, J approximately corresponds to the typical amplitude of the 
Josephson current expected in a single sample. The relation 
(J) ~ has different meaning for SFS and S/HM/S cases. In 
SFS junctions, (J) = at the transition points is the result 
of the ensemble averaging and the Josephson current remains 
finite in a single sample. The characteristic temperature and 
length of a ferromagnet at the — vr transitions vary from one 
sample to another. In S/HM/S junctions at Vex — 2.5t, how- 
ever, (J) = means vanishing Josephson current even in a 
single sample'^ because (J) ^ 5J = 0. 

The origin of large fluctuations of the Josephson current 
can be understood by considering the behavior of the pairing 
function in a ferromagnet. The pairing function in Eq. (|3} can 
be decomposed into four components. 



^ W 3 



(5) 



where /o(/3) is the pairing function of the spin-singlet (spin- 
triplet) pairs with the spin structure of (|ti) — (+) UT))/'s/2, 
respectively, and the pairing function of Ifj) (lii)) is given 
by /tT = - fi (III = */2 + /i)- In Figs- fib) and (c), 
we show (/o) and 6fo as a function of position j in a dif- 
fusive ferromagnet, where /b is the pairing function in bulk 
superconductor, uin is fixed at 0.02Ao, Vs — 0, and ip — 0. 
The junction interface and the center of a ferromagnet corre- 
spond to j = 1 and j — 37, respectively. In SNS junctions 
(Fig.l^'), (/o) is larger than 6fo and very weakly decays with 
j. The spin-singlet Cooper pairs exist everywhere in a normal 
metal. Near the interface, 6fo is slightly suppressed due to 
the tight contact to the superconductor. On the other hand, in 
SFS junctions in (Fig. the average (/o) decreases expo- 
nentially with j according to exp{—j /£^h) as indicated by a 
broken line. The fact that 6fo remains finite at the center of 
a ferromagnet means that the spin-singlet pairs penetrate far 
beyond even though (/o) ^ there. 

In Fig.|2ld), we show the pairing function for three differ- 
ent realization of disorder in SFS junctions. The pairing func- 
tions are in phase near the interface (j < ^h), whereas they 
are out of phase far from the interface. We obtain the rela- 
tion (5 A) oc e~^/^^ with = ^Djlu^n in agreement with 
Ref. Y% Thus we conclude that spin-singlet Cooper pairs do 
exist in a single sample of ferromagnet even for j ^ ^/j and 
the mesoscopic fluctuations of the pairing function provide the 
origin of the large fluctuations in the Josephson current. In 
S/HM/S junctions for = 2.5t as shown in Fig.|2jb), both 
(/o) and 6 fa vanish for j » 1, which indicates the absence of 
spin-singlet Cooper pairs in a half metal. 

The relation ( J) <C SJ is the characteristic feature of the 
Josephson current in diffusive SFS junctions. This feature, 
however, is drastically changed by the spin-flip scattering at 
the interfaces. In Figs.|5](a) and (b) we show (J) and 6J vs 
Vs for Vex/t = 1 and 2.5, respectively. In both cases (a) and 
(b), we find that |(J)| > 6J for Vs > OM. The Josephson 
current becomes self-averaging in the presence of the spin- 
flip scattering. The reason can be explained by calculating the 
pairing functions of equal-spin pairs shown in Figs. |3jc) and 
(d), where fi, is plotted as a function of position j. Here we 
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FIG. 3: (Color online) (a) The Josephson current and its fluctuations 
at T = O.lTc and ip — n/2 as a function of the interface spin-flip 
scattering Vs for Vex/t — 1 and (b) for Vex/t — 2.5. The pairing 
functions versus position j in a ferromagnet (c) and in a half metal 
(d) at Vs = OAt, V3 = and tJ„ = 0.02Ao. In (d), differences 
among (fn), Sfo and S fs at large j become small. 



show Sfo and Sf^ instead of (/o) and (/s) because the ensem- 
ble averages are much smaller than their fluctuations. The fast 
decay of Sfo and Sf^ is determined by strong spin polariza- 
tion. In both cases (c) and (d), becomes larger than Sfo 
and 5/3 because the pairing function /n does not change sign 
for various impurity configurations. Thus the averaged quanti- 
ties become larger than their fluctuations. Thus the Josephson 
cuiTent becomes self-averaging as shown in Figs.lSja) and (b). 

Finally we address an unusual symmetry property of the 
Josephson current in S/HM/S junctions. In Fig.|3(a), we show 
(/ll) as a function of w„, where j = 37, Vs = 0.2t, (p = 0, 
and Vex = 2.5t. For comparison, we also show (/o) on the 
normal side of a SNS junction. The pairing function (/□) in 
a normal metal is the even function of cj„, whereas (/||) in 
a half metal is the odd function of ojn Ref. 6. The pairing 
function obeys the Pauli's rule 



/„„(r,r') = -/-.„(r',r), 



(6) 



where / denotes the transpose of / meaning the interchange of 
spins. It is well known that ordinary even-frequency pairs are 
classified into two symmetry classes: the spin-singlet even- 
parity and the spin-triplet odd-parity one. In the former case, 
the negative sign arises due to the interchange of spins, while 
in the latter case due to r ^ r' . In the present calcula- 
tion, all components on the right hand side of Eq. Q have 
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FIG. 4: (Color online) (a): Dependences of the pairing functions 
on u!n- (b): The local density of states at j = 37 in a half metal at 
Vs = 0.5t and in a normal metal at Vs = 0. 



the s-wave symmetry. The pairing functions are isotropic in 
both the real and momentum spaces due to diffusive impu- 
rity scattering. As a result, /|| must be the odd function of 
u!n to obey the Pauli's rule. Both even- and odd-frequency 
pairs are mixed in ferromagnets as shown in Fig. |3c). The 
fraction of odd-frequency pairs depends on parameters such 
as the exchange potential and the spin-flip scattering. On the 
other hand, in a diffusive half metal all Cooper pairs have the 
odd-frequency character, which causes drastic change in the 
quasiparticle density of states. 

The density of states is given by N{E,j) — 
-iwI]m=iImTrG'£;+i^(r,r), where 7 is a small imagi- 
nary part chosen to be 0.05 Aq in the following. In Fig.|4jb), 
the local density of states (LDOS) at j — 37 is shown, 
where ip = and A^o is the density of states in the normal 
state at Vex = 0. For comparison, we show LDOS on the 
normal side of SNS junction with Vs = which has a 
minigap at E < Exh ~ O.SAq, where Eth is the Thouless 



energy. In contrast to that, LDOS in a half metal has a 
peak at the Fermi energy, it's width is characterized by 
Exh- This peak is generated at t he sp in active interface by 
the mechanism discussed in Ref 20 and is tranfered into 
a half metal due to long range property of odd-frequency 
spin-triplet even-parity pairing function. The peak is much 
stronger than the enhancement of the LDOS found in weak 
ferromagnets^iiSi^iii^. In addition, in a half metal the peak 
shape is almost independent of position, while in the SF 
junctions^ the LDOS has an oscillatory peak/dip structure 
at E = which rapidly decays with the distance from the 
SF interface. Therefore the large peak at = in LDOS 
is a robust and direct evidence of the odd-frequency pairing 
in half metals. To test the existence of such peculiar pairing 
state, the scanning tunneUng spectroscopy could be used. 

In conclusion, we have studied Josephson effect in super- 
conductor / diffusive ferromagnet / superconductor junctions 
by using the recursive Green function method. The Joseph- 
son cuiTent in these junctions basically is not self-averaging 
because the spin-singlet Cooper pairs penetrating into feiTO- 
magnets far beyond cause the large fluctuations of the pair- 
ing function. In the presence of the spin-flip scattering at the 
interfaces, the equal-spin odd-frequency pairs drastically sup- 
press the fluctuations. When ferromagnets are half-metallic, 
all Cooper pairs have the odd-frequency property. As a result, 
the low energy peak in the quasiparticle density of states in a 
half metal exists at the distances far beyond from the inter- 
face and could be probed by scanning tunneling spectroscopy. 
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